Calculations of electric fields for radio detection of Ultra-High Energy particles 



Daniel Garcia- Fernandez, Jaime Alvarez-Mufiiz, and Washington R. Carvalho Jr. 
Depto. de Fisica de Particulas & Instituto Galego de Fisica de Alias Enerxias, 
Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain 

Andres Romero- Wolf 
Jet Propulsion Laboratory, California Institute of Technology, 
4800 Oak Grove Drive, Pasadena, Caltforma 91109, USA 



(N 

O 

(N 
> 

o 

O: 
6 



> 

o 
d 



X: 



Enrique Zas 

Depto. de Fisica de Particulas & Instituto Galego de Fisica de Altas Enerxias, 
Universidad de Santiago de Gompostela, 15782 Santiago de Compostela, Spain 

The detection of electromagnetic pulses from high energy showers is used as a means to search 
for Ultra-High Energy cosmic ray and neutrino interactions. An approximate formula has been 
obtained to numerically evaluate the radio pulse emitted by a charged particle that instantaneously 
accelerates, moves at constant speed along a straight track and halts again instantaneously. The 
approximate solution is applied to the particle track after dividing it in smaller subintervals. The 
resulting algorithm (often referred to as the ZHS algorithm) is also the basis for most of the sim- 
ulations of the electric field produced in high energy showers in dense media. In this work, the 
electromagnetic pulses as predicted with the ZHS algorithm are compared to those obtained with 
an exact solution of the electric field produced by a charged particle track. The precise conditions 
that must apply for the algorithm to be valid are discussed and its accuracy is addressed. This 
comparison is also made for electromagnetic showers in dense media. The ZHS algorithm is shown 
to describe Cherenkov radiation and to be valid for most situations of interest concerning detec- 
tors searching for Ultra-High Energy neutrinos. The results of this work are also relevant for the 
simulation of pulses emitted from air showers. 



I. INTRODUCTION 

The study of Ultra High Energy Cosmic Rays 
(UHECR) and neutrinos (UHEz^s) is currently a high pri- 
ority in Astroparticle Physics with many experimental 
efforts being dedicated to these two related areas of re- 
search. UHECRs are routinely detected through the Ex- 
tensive Atmospheric Showers (EAS) they produce when 
interacting in the atmosphere. Despite the recent ad- 
vances in the measurement of the flux of UHECRs P, Q , 
their primary composition remains unknown |3| , and 
this is one of the main obstacles to extract precise con- 
clusions on their origin. There are strong reasons to be- 
lieve that UHEi^s should be produced in the interactions 
of UHECRs with the material surrounding the sources, 
and/or in their propagation through the observed Cosmic 
Microwave Background radiation Q. However their de- 
tection has not yet been achieved. Efforts are being made 
to improve the experimental situation in both fields. The 
challenge is to instrument sufficiently large target vol- 
umes to compensate for the low cross section and low 
fluxes for UHE:^ detection, and to find measurements 
that provide large aperture at the highest energies and 
help to constrain the primary composition of UHECRs. 
The radio technique is being explored in both fields. 

As early as in 1962, G. Askaryan proposed to detect 
UHECRs and UHEi^s by observing the coherent radio 
pulse from the excess of electrons in a shower developing 
in a dense, dielectric and nonabsorptive to radiowaves 
medium Q. Soon after, pulses were observed in coin- 



cidence with air shower arrays 0, 01 • The emission is 
coherent in wavelengths which are large compared to the 
characteristic size of the electric charge and current dis- 
tributions associated to the induced showers. The tech- 
nique has been receiving a lot of attention in the field of 
Astroparticle Physics in the last decade because of the 
relatively low cost of the antennas needed for the detec- 
tion systems. Also, coherence implies that the pulse en- 
ergy scales with the square of the primary energy which 
favours long range detection, a requirement to achieve 
large areas and volumes both for UHECR and UHEz^ de- 
tection. Indeed, quadratic scaling has been confirmed 
in accelerator experiments [9-12] leading to very strong 
pulses associated to UHE energy showers. 

A large number of initiatives have been made or are 
currently in development or planning stages. In dense 
media neutrino-induced showers are less than a meter in 
width and full coherence is expected to be mantained up 
to the GHz range (see for instance [ill). ^ variety of past 
and present experiments search for pulses at those fre- 
quencies produced by neutrinos in ice the moon 
regolith [16l - l21] or salt domes [11| . The largest and most 
promising experiments are in planning stages looking at 
ice [25, ^235 . Electromagnetic pulses emitted in the MHz- 
GHz frequency range by EAS induced by UHECR are be- 
ing measured in the hope of using this complementary in- 
formation to constrain its composition and/or to develop 
new cost effective detection systems 24, 25). In addition 
the ANITA balloon flown antennas, initially devised to 
search for neutrinos, has recorded coherent pulses up to 
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the GHz range [2^ consistent with EAS induced by UHE- 
CRs and there are plans for future developments (273 . 

The success of the radio technique requires an accu- 
rate and computationally efficient calculation of the ra- 
dio emission properties of UHE showers. It is necessary 
to perform this in an efficient way since the number of 
particles in a shower at EeV energies is > 10^. Sim- 
ulation techniques for evaluating pulses in dense media 
have been used for more than 20 years [l^ [28l - l4]| . The 
shower is simulated to obtain, for every particle track, the 
information needed to calculate its contribution to the 
electric field. The formula used for this purpose stems 
from an approximate solution of Maxwell's equations in 
the Fraunhofer limit [8] adapted for simulation purposes 
in The electric field due to a short charged particle 
track, assumed to be travelling at a constant speed, can 
be approximated by two terms which correspond to the 
start and end of the track. The resulting electric field 
depends on the particle speed and the angle between the 
line of sight from the track to the observer [l^, ■ To 
calculate the emission in a shower, tracks are chosen so 
that the particle velocity can be approximated to be con- 
stant, and all track contributions are added taking into 
account interference effects. 

Although the approximate formula is sufficient for 
many practical applications, its range of validity is lim- 
ited in frequency and position of the observer with re- 
spect to the track. It is possible to extend its range of va- 
lidity by subdividing each track in sub-intervals. The re- 
sulting algorithm (often referred to as "ZHS algorithm" ) 
is easy to implement and fast enough for the simulation 
of particle showers. This is convenient since it has al- 
lowed the simulation of pulses in the Fresnel region for 
neutrino detection [s^l and in measurements of EAS [13] . 
Although the range of applicability of the algorithm is en- 
hanced when used in this way, it is not obvious that the 
sum of the subcontributions correclty accounts for all the 
radiation in regions close to the emission source. 

The object of this article is to study the conditions for 
the ZHS algorithm to be valid and to establish its accu- 
racy. For this purpose we obtain in Section II an exact 
solution to the problem of a charged particle instanta- 
neously accelerating to a constant speed and stopping 
abruptly after a discrete time interval [i^. This allows 
us to establish the precise conditions necessary to turn 
this solution into the basic formula of the ZHS algorithm. 
In Section HI we compare the exact solutions for an infi- 
nite and a finite track to the result of the ZHS algorithm, 
and stress the compatible interpretation of the radiation 
regime in terms of Cherenkov radiation for both cases. 
In Section IV we compare the exact solution of the sin- 
gle track problem in nearby regions with the results of 
applying the ZHS algorithm with track subdivisions to 
test its validity and accuracy. Section V is devoted to 
discussing the ZHS algorithm in relation to other ap- 
proximations made to calculate the pulses emitted from 
showers. Section VI presents the summary and conclu- 
sions of our work. 



II. ELECTROMAGNETIC FIELD OF A SINGLE 
CHARGED PARTICLE TRACK 

Let us assume an electron ejected from an atom at time 
t = ti that travels at constant speed v through a medium 
along a finite track until it is absorbed by another atom 
at time i2- Neglecting the movement of the atoms, we 
can model the electric current associated to the electron 
as, 

J(x,t) = -ev(5(3)(x-xo - vt) e(t-ii)e(i2 -t) (1) 

where e = |e| is the charge of a positron, x(t) is its po- 
sition and Xq an arbitrary reference position. The step 
0-functions account for the fact that the electron only 
moves in the time interval (ii, ^2)- 

In a dielectric medium with permitivity e and magnetic 
susceptibility /z. Maxwell's equations for the vector po- 
tential in the frequency domain A(x, w) can be written 
as B]: 

V^A(x, tj) + ^ea;^A(x,a;)- 

V[V • A(x, w) — iefiLU 0(x, w)] = — /iJ(x, w), 

where (/)(x, w) is the Fourier transform of the scalar 
potential and we use the following convention for 
the Fourier-transform of a function f{t): /(w) = 
J^^dt e*"*/(0- principle e and /x can depend on 
frequency and our results below would be equally valid, 
but we drop the explicit dependence of e and /i with uj 
for simplicity. 

We use the Lorenz gauge condition which implies that 
Eq. ^ for the vector potential becomes: 

W^A + k^A = -fiJ, (3) 

with k ~ nuj / c and n the refractive index of the medium. 
The Lorenz gauge condition in the frequency domain j44| : 

V • A(x, oj) = ie/iw (/)(x, w), (4) 

implies that the scalar potential for non-zero frequencies 
is entirely determined by the divergence of the vector and 
only the vector potential is needed to calculate the field. 

A. Exact solution 

The solution of the Helmholtz equation in Eq. (jS]) is 
standard physics and can be obtained, using Green's 
method, as an integral over source positions [4J]: 

A(x,a.) = d=^x'^— -^J(x',..) (5) 

where from now on, x will denote the position of the 
observer and x' the position of the source. 

For simplicity we rewrite the current in Eq. ([T]) as: 

J(x,t) = qv Z{t) P(x,t) z (6) 
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with q = -e; Z{t) = Q{t - ti)Q{t2 - t) and P(x,t) = 
,5(3) — vtz — zqz), where we assume without loss of gen- 
erality that the electron travels along the z-axis which is 
parallel to z. The equations that follow below are clearly 
valid for any continuous and differentiable functions Z(t) 
and P(x, t). They are also valid for the case considered 
in Eq. ([I} which can be obtained as a limiting case of 
suitable continuous and differentiable functions. 

Transforming this definition of the current to the fre- 
quency domain and substituting into Eq. ([5]) for the vec- 
tor potential we obtain: 



An J |x — x'l 

(7) 

In the Lorenz gauge, the only non-zero component of 
the vector potential for a charge moving along z is the 
z-component. The scalar potential (j) is obtained through 
Eq. (|4]). For this purpose we need to obtain V • A: 



V • A(x, w) 



dz Att 

ifc|x- 



qv 



J d^x'di' Z{t')P{jc',t') 
{z-z') 



ik 



The electric field E(x, w) can be obtained from the 
scalar and vector potentials 44] as: 



E(x, uj) = — V(/)(x, uj) + iwA(x, uj) 



(8) 



Due to the cylindrical symmetry of a charged particle 
track we calculate only the radial (Ep) and z {E^) com- 
ponents of the field. Since A only has z-component, then 
the radial component of the field is given by: 



E,{^,Lo) = I — / At' e*"* ^ X 

Lo Attc 

2 {Z- Zq- vt'Y (Z- Zq- vt'Y 



(12) 



iuj J^qv / At' e'^' — 
An r 

where r = r{t') and h — b{t') are both functions of the 
source time t' which is defined as, 



and 



-{t') = |x - x'l = ^/p^ + (z - zo - vt 



h(t') =ik^ 



l\2 



r{t') 



(13) 



(14) 



Eqs. PT|) and provide an exact solution for the 
electric field of a finite track. In general they do not have 
an analytical form and a numerical integration needs to 
be performed to obtain the field. For the calculations 
in this paper, we have divided the integration interval 
and applied Simpson's rule, increasing the number of di- 
visions until the integral converged. Under certain con- 
ditions, however, we can give analytical approximations 
for relevant physical situations. 

In the following we show that the basic expression used 
in the ZHS algorithm is a particular case of Eqs. (fTTj) and 
(|12p under certain approximations. 



Ep{yi,u) = 



^0(x,cj) 
dp 



V • A(x,w) 



(9) 



where p — y^zJ^ + j/^ is the radial coordinate, and dp 
denotes d/dp. 

The z— component of the field is given by: 

E;,{-K,uj) = ~dz(l){:x.,Lo) + ia;A^(x,iu) 

„ V • A(x, uj) . , X „N 

iujpe 

with dz denoting d/dz. 

Performing the derivatives in Eqs. © and pH)) we ob- 
tain: 



qv 1 



Ep{^,Lo) ^i^i^ I dt' e*"* 

U! Ane Ji^ r-^ 



p X (z — Zq — vt') X 



b b 



(11) 



B. The ZHS expression 

A simple expression for the approximate calculation 
of the electric field from a single charged particle track 
moving at constant speed was found in 13]. In this sec- 
tion we derive the expression used for the ZHS algorithm 
from the exact solution and compare the electric field 
as obtained in both the exact calculation and with the 
ZHS expression. This allows us to establish under which 
circumstances the formula gives a good account of the 
electric field. 

The ZHS algorithm can be obtained from Eqs. pT|) 
and (|12p if the following set of conditions are fulfilled: 



1. The observer is in the "far field" zone i.e. 

kr > 1 (15) 

2. The Fraunhofer approximation holds. This can be 
stated as a condition for the phase factor to be 
approximated as: 



and, 



kr = fc|x — x'l w k[R — v{t — to) cos ( 



(16) 
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where R is the distance from the observation point 
to a reference point along the track where the par- 
ticle is located at a reference time ^o, and 9 is the 
angle between the particle track and the direction 
from the reference point to the observer. This ap- 
proximation holds provided the parameter rj ^ 1 
with 77 defined as: 



k[v{t - to)] 
R 



(17) 



This condition should be fulfilled at any time t from 
ti to ^2- A more commonly used and nearly equiv- 
alent form of this condition WSl is: 



kl? 



■ sm 



« 1, 



(18) 



where L = w(t2 — ii) is the length of the track. This 
condition is necessary to ensure that the second 
and higher order terms for the phases i[ujt + kr) 
in Eqs. pT|) and ((T^ have no relevance, even when 
the sum of the leading and first order terms in the 
Taylor expansion of the phases is zero, as it occurs 
for observation at the Cherenkov angle 49] defined 
as cos6'c; = 1/Pn with fi = v/c. 

3. Finally the distance to the observer appearing in 
the denominators of several terms of Eqs. PT|) and 
([T^ must be approximated as: 



1^1 
HT)^ R 



(19) 



over the length L of the track, where R is the dis- 
tance to a reference point along the track (in the 
ZHS algorithm the mid-point of the track is se- 
lected) . The error when making this approximation 
is of order L/R. 

In particular, the condition in Eq. (llSp implies: 



b{t') « ik 



b--Kb 



b^^ — Kh' (20) 



With these approximations the radial component of 
the field in Eq. (fTTj) becomes: 



£; ~ J_ ^ Ukf sin 6* cos 6* e*"^-^*" 
LJ Aire R ^ ' 



dt' e'("-'"^^*' 



(21) 



where we have used: 



r ^ R 



sin 9, 



and 



z — Zq — vt Z — Zq — vto 
r ^ R 



cos9. 



(22) 



(23) 



Eq. ([21]) can be easily integrated yielding: 

,. ^ikR 



Ep = —iq u! - — w sin cos 



47r 

gi(w-k-v)t2 _ g 



R 

i{uj—'k.-v)ti 



z(C(J — k • v) 



(24) 



If we make to = ti this becomes the expression for the 
radial field as used in the ZHS algorithm [l3| except for 
a factor 2 due to the Fourier transform convention used 
in [13]. 

Similarly applying the approximations in Eqs. ([16]), 
([T9l) and ([20|) to Eq. ([12]) for the z-component of the 
field, and using that kR ^ 1 fc^ ^ k/R, it is straight- 
forward to show that: 



-iLo^ - — e*" ^ cos^ I 



uj Att R 
47r i? 



di e'^""""-^)* 



ti 



(25) 



which can be cast as: 



/i 2 
Ez = iq LJ — V sm I 



„ikR 



R 



/ di e*("-''-^)* (26) 



Performing the integral, and taking to — ti, the ZHS 
formula is recovered: 



E, 



iq uj 



— V SIT? 9 



JkR 



R 



gi(w-k-v)t2 „ gi(cj-k-v)ti 
l(aj — k • v) 



(27) 



C. The ZHS algorithm 

To calculate the electric field of the pulse emitted from 
a current distribution, such as that produced in a high 
energy shower, the ZHS algorithm uses the ZHS expres- 
sions in Eqs. ([24)) and ([27[) to calculate the emission from 
all the charged particle tracks. The final result is ob- 
tained adding up all the contributions. The value of R 
used for the phase factor in each particle track is the dis- 
tance from the first point of the track to the observation 
point. This definition is consistent with the convention 
to account for the phase change between emission arising 
from the start and end points of the track. The actual 
value of R used for the denominator is the distance be- 
tween the midpoint of the track and the observer which 
is for all practical purposes the same value when the ap- 
proximation is valid. The algorithm used in alternative 
simulation programs is similar but different in these tech- 
nical details [45 1. 

Naturally it is possible to divide any charged particle 
track into arbitrarily small subtracks in order to make 
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the computation more accurate. This will extend the 
range of validity of the approximation. It is interesting to 
note that it makes absolutely no difference to subdivide 
the track of a uniformly moving charge when observed 
in the Fraunhofcr limit because the term associated to 
the end of one subtrack cancels the term associated to 
the beginning of the next subtrack. This is because in 
this limit the differences in R between adjacent subtracks 
are arbitrarily small. However in practical situations this 
cancelation is not exact if R is allowed to change from a 
track to the next one according to geometry. In the origi- 
nal ZHS program @ the tracks were not subdivided but 
soon it was realized that a more accurate result was ob- 
tained by subdividing all tracks at every point there was 
a discrete interaction in the simulation program ^35, '48'] . 
This was the main modification that was ever made to 
the original ZHS Monte Carlo [l^, and has been effec- 
tive since then. With this subdivision the distribution of 
tracks for a shower in ice has a peak at about L ~ 1 mm. 



the speed of light, then ^(a;) is imaginary and the particle 
radiates. The latter case corresponds to pure Cherenkov 
radiation [49^ . 

With the help of the asymptotic forms for the Bcssel 
functions, we can obtain the limits of Eqs. (|28p and ((29)) 
when \up\ <^ 1 i.e. for small distances to the track com- 
pared to the radiation wavelength. Conversely we can 
also obtain the limits when \up\ 3> 1 i.e. for large dis- 
tances compared to the wavelength. If \up\ ^ I the Ki 
Bessel function dominates over the Kq and only the ra- 
dial component of the field Ep matters. In this case, 



lim \E„ 



\q\ 1 



2TTev p 



(32) 



and a 1/p dependence with distances is obtained, as well 
as no dependence with frequency. 

If I up I ^ 1 and keeping in mind that when v > c/n 
the argument u is imaginary the fields can be written as: 



III. CHERENKOV RADIATION 



lim \Ep\ = 

Lip— >-±ioo 



2TTev 



V ""P 



(33) 



The exact electric field for a charged particle track 
derived in Eqs. (ITTt and must account for every 
single feature of the electric field, since no approxima- 
tions have been made. In particular, it must reproduce 
Cherenkov radiation which is the only radiation emitted 
by a charged particle moving at constant speed, when 
V > c/n in the limit of an infinite track. An analytical 
solution for the electric field produced by such particle 
can be obtained and it is given in [i^ (chapter 4). The 
radial p and z— components of the fields given in (49l | 
converted to SI units, and using the Fourier transform 
convention adopted in this work, can be written as: 



Ep{p,z,uj) 



2TTev 



^''uKi{up) 



(28) 



lujfiq 

f 

27r 



Koiup) , (29) 



where p is the radial coordinate, Kq and Ki are the mod- 
ified Bessel functions of the second kind, and u = u{lo) 
is a function that can take two different values depend- 
ing on the magnitude of the particle speed, v < c/n or 
V > c/n (subluminal or superluminal regime). In our 
convention: 



c w 

U < — U(UJ) = — 

n V 



V > — ^ u(lu) 
n 



-i — 

V 



1 



(30) 



(31) 



If the particle travels below the speed of light in the 
medium, the argument uicd) of the Bessel functions is 
real and the particle does not radiate as shown in p9l |. 
On the contrary, if the speed of the particle is larger than 



lim \E^ 

up-^:tioo 




— (34) 



In this case the field is proportional to y/uj/p. This is in 
agreement with (46j where the same behavior is deduced 
using simple arguments of energy conservation through 
a cylindrical surface surrounding the track. 

In Fig. [T] the Fourier components of the modulus of 
the electric field for an infinite track as obtained from 
Eqs. (pS)) and (|29p are shown, for different frequencies, 
as a function of p, the radial distance to the track. The 
particle speed is u ~ c > c/n travelling in homogeneous 
ice with refractive index n — 1.78. Under these circum- 
stances the quantity \up\ can be approximated as: 



\up\ 



of " ] 

V 100 MHz/ Vim/ 



(35) 



At large distances to the track when \up\ ^ 1 the fields 
shown in Fig. [1] scale with distance as l/^/p and with 
frequency as ^/uj, in agreement with the aysmptotic field 
components in Eqs. I\'3'3\i and ([M)) . This behavior takes 
place when p > 0.1, 1 and 10 m for frequencies 1^=1 
GHz, 100 MHz and 10 MHz respectively, in agreement 
with Eq. (1551) . as can be seen in Fig. [TJ 

As the distance to the track decreases and the con- 
dition \up\ <C 1 starts to be valid, the field behaves as 
l/p and does not depend on frequency as expected from 
Eq. (|32|) . As can be clearly seen in Fig. [1] the transition 
from the l/p behaviour to \fijjjp occurs at a distance 
that depends on frequency because \up\ involves the fre- 
quency (Eq. 1551) . For instance at distances p < 0.01 m 
the condition \up\ <C 1 applies for both i/ = 100 and 10 
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MHz. The Fourier component of the field scales with \ j p 
and has the same value for the two frequencies as seen in 
Fig. [U while this is not the case for a frequency of = 1 
GHz. 




0.1 1 10 100 1000 

P [ml 



FIG. 1. Fourier components of the electric field modulus 
as a function of distance to the particle track for an infinite 
track as obtained from Eqs. (|28|) and (|29|) (solid line), and 
for a track of length L — f 200 m as obtained with the exact 
formulas derived in this work (Eqs. (1111) and (112^ 1 (open cir- 
cles) and with the ZHS algorithm (Eqs. ((24)1 and ((271)) (open 
squares). From top to bottom, the observation frequencies are 
1 GHz, 100 MHz and 10 MHz. The 1/p and l/^p regimes 
are apparent. 




V [MHz] 



FIG. 2. Top panel: Fourier components of the electric field 
modulus as a function of frequency for an infinite track as 
obtained from Eqs. (|28p and (|29p (solid line), and for a track 
of length L = 1200 m as calculated in this work from Eqs. pip 
and (|12p (open circles). The observer is placed at a lateral 
distance to the infinite track p — 8.27 m. Also shown is the 
modulus of the field for the same finite track (L — 1.2 km) as 
obtained with the ZHS algorithm, Eqs. (pi)) and ^} (open 
squares). At high frequencies, the field behaves with ^/uu (see 
text for explanations). Bottom panel: Relative difference (in 
%) between the solution for an infinite track and the exact 
solution for a finite long track (open circles) and between the 
solution for an infinite track and that obtained with the ZHS 
algorithm (open squares). 



In Fig. [2] the modulus of the field for a charged parti- 
cle in an inhnitc track - as obtained from Eqs. ((28)) and 
P9[) - is shown as a function of frequency for an observer 
at a fixed radial distance. At large enough frequencies 
so that the condition \up\ ^ 1 applies, the field scales 
as -y/w as expected from Eqs. ((33|) and ([M)) . while it is 
constant with frequency for small enough frequencies so 
that \up\ ^ 1 as predicted from the asymptotic Eq. ((5^. 
More quantitatively, since the observer in Fig. (2] is lo- 
cated at /9 ~ 10 m, the field should behave as y/uj for 

Jc; 10 MHz (applying Eq. ((55)) ). This is approximately 
the case as can be seen in Fig. (21 

The result of the exact calculation for a track of length 
L = 1.2 km is also shown in Figs. [Tj and [2j The agreement 
between both calculations is excellent because the finite 
track is long compared to observation distance. This 
confirms that Eqs. ((TT)) and ((12)) must also account for 
Cherenkov radiation. As can be appreciated in Fig. (21 the 
exact results for finite and infinite tracks differ for wave- 
lengths larger than the length of the track - frequencies 
typically below vq ~ {c/n)/X with \ L = 1.2 which 
gives V <0.l MHz. 



IV. COMPARISON OF THE EXACT 
CALCULATION AND THE ZHS ALGORITHM 

We have numerically evaluated the exact expressions 
for the z and p components of the electric field in 
Eqs. ()lip and ((T^ at different frequencies and observer 
distances, for a single tracks of different lenghts, and for 
the tracks constituting a shower in ice as obtained in full 
simulations performed with the ZHS Monte Carlo code 
p^ . In this section we present the results of this com- 
parison. 



A. Fourier components of the electric field for a 
single track 

The applicability of the ZHS expressions in Eqs. ((24)) 
and ((?f)) relies on the conditions 1 — 3 in Section IIIBI 
Eqs. ((T7)) and ((TO)) are easily fulfilled for any wavenumber 
k simply by dividing the particle track in a sufficiently 
large number of subtracks. Once this is guaranteed, the 
ZHS expression is applied to every sub-track and the elec- 
tric field is obtained adding the corresponding contribu- 
tions. The validity of this procedure will be numerically 
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confirmed below when comparing the exact calculation 
with the electric field obtained using the ZHS algorithm 
in the manner just described. However the condition 
kr 1, Eq. dTSl) . does not depend on the size of the 
track and cannot be enforced by applying the procedure 
outlined above. As a consequence fcr ^ 1 is an intrinsic 
limit to the range of observing frequencies and distances 
in which the ZHS algorithm gives accurate results as will 
be shown in the following. 




V [MHzJ 

FIG. 3. Top panel: Fourier components of the electric field 
modulus for a single particle track as obtained with the exact 
calculation Eqs. and (|12|l (solid lines) and with the ZHS 
algorithm Eqs. (f24)l and (f27|l (open symbols). The length 
of the track is L = 1.2 10""^ m and the field is shown for 
observers at distances (from top to bottom lines) 7? = 0.1 m, 
1 m, 10 m and 100 m, with respect to the center of the track, 
and placed at the Cherenkov angle. Bottom panel: Relative 
difference (in %) between the exact solution and that obtained 
with the ZHS algorithm for the same distances. 



In Fig. [3] we compare the Fourier components of the 
electric field modulus for a single particle track as ob- 
tained with the exact calculation, Eqs. (fTT|) and ([T2l) . to 
that obtained with the ZHS algorithm. The length of 
the track is chosen to be small L = 1.2 10~^m (close 
to the peak value of the distribution of track lengths in 
the standard ZHS code) . To test the validity of the ZHS 
algorithm, we have calculated the spectra for observers 
at different distances (i?) measured with respect to the 
center of the track and placed at the Cherenkov angle. 

The condition kR ^ 1 in ice with refractive index 
n — 1.78 can be cast as: 



( - ) 

V 100 MHz/ Vim 



> 1 



(36) 



For observers at distances R = 100, 10, 1 and 0.1 m 
from the particle track, the condition in Eq. (j36l) is ful- 
filled as long as > 1, 10, 100 MHz and 1 GHz respec- 
tively. The ZHS algorithm is expected to reproduce the 



results of the exact calculation in this range. This can be 
seen in Fig. [31 The relative difference between the ZHS 
and exact calculations is less than ~ 2% at the fron- 
tier of the validity range in the explored frequency and 
distance space. The accuracy can be however orders of 
magnitude better. If kR > 37 is enforced for instance 
the corresponding relative difference is below ^ 0.01%. 
This condition is satisfied for i? > 10 m and i' ^ 100 
MHz what corresponds to a range of frequencies and dis- 
tances typically encountered in experiments that search 
for neutrino induced radio transients. 

What is striking is that these conclusions also apply 
in the "near field" provided the track is subdivided in 
sufficiently small sub-tracks. This can be seen in Fig. |4] 
comparing the exact and ZHS results for a L = 1.2 m 
track. The same range of validity is obtained because 
conditions ((T6| , (flT)) and (fT9)) are guaranteed by reducing 
the length of the sub-tracks. 



10" 




Exact calculation 

100 1000 



V [MHz] 



FIG. 4. Same as top panel of Fig.[3]for a single particle track 
of length L — 1.2 m. 



For fixed frequencies the condition in Eq. (|36|) only ap- 
plies at sufficiently large distances R to the track. This 
is illustrated in Fig. [5] for a track of length L = 1.2 m. 
The Fourier components at frequencies of ly — 10, 100 
MHz and 1 GHz are in agreement with the exact calcu- 
lation respectively at i? > 10, 1 and 0.1 m as expected. 
Since the typical distance between antennas in experi- 
ments such as the Askaryan Radio Array (ARA) [22] is 
~ 10 — 100 m, we expect the results of the ZHS algo- 
rithm obtained through the procedure outlined above, to 
be accurate enough in most practical situations. 

It has been questioned whether the ZHS algorithm re- 
produces Cherenkov radiation from a single charged par- 
ticle track 45|. In Figs. [T] and [5] the algorithm is shown to 
be in very good agreement with the exact solution for a 
1.2 km track as long as kR ^ 1 as explained above. The 
emission from such track is in turn practically equivalent 
to the Cherenkov emission from an infinite track (also 
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displayed) and it follows that the ZHS calculation must 
account for Cherenkov radiation. 




0.01 0.1 1 10 100 1000 

Rim] 



FIG. 5. Top panel: Fourier components of the electric field 
modulus aX u ^ I GHz, 100 MHz and 10 MHz (from top 
to bottom lines) for a single particle track as a function of 
distance to the track, as obtained with the exact calculation 
Eqs. (fTT|) and (fT2l) (solid lines) and with the ZHS algorithm 
Eqs. ((211 a-nd (f27|) (open symbols). The length of the track 
is L = 1.2 m and the observers are placed at the Cherenkov 
angle. Bottom panel: Relative difference (in %) between the 
exact solution and that obtained with the ZHS algorithm for 
the same frequencies. 



1. Behaviour of the field with frequency 

As can be seen in Fig. [3] (solid lines) the exact solu- 
tion for the modulus of the electric field scales linearly 
with frequency provided that kR ^ 1. As a result the 
electric field can be approximated with the ZHS expres- 
sions, Eqs. (j24p and (p7)) . For an observer close to the 
Cherenkov angle as in Fig. [3l the factor (a; — k • v) <C 1 
and the term in brackets in Eqs. (|24p and (f27| can be 
approximated by: 

"gi(i^-k-v)t2 _ gi(cj-k-v)ii " 



1 + iiuj — k • v)i2 — 1 — «(w — k • v)ti 



making the dependence of the field with lo apparent. 
Clearly the field cannot grow indefinetely with frequency. 
This apparent "ultraviolet divergence" is only an arti- 
fact of considering the unrealistic medium in which the 
permittivity e(w) is constant with frequency. In a physi- 
cal medium absorption at high frequencies will tame the 
growth with frequency of the electric field. 

When kR < 1 the field behaves with frequency as uj~^ 
as can be also seen in Fig. |31 In the model of a charged 



particle at rest for t <ti, moving with a speed v between 
t = ti and t = t2, and becoming again at rest for t>t2, 
the Coulomb field dominates at small distances to the 
track and/or low frequencies. The field can be modeled 
as: 

j5i cx ^ Q{h,obs - t) (38) 

for t <ti and 

£^2 OC ^ e(i - t2,obs) (39) 
^2 

for t > t2, where ri and r2 are respectively the distances 
from the charge to the observer at times ti and t2, and 
ti^obs and t2,obs denote the instants of time at which the 
Coulomb field arrives at the observer. The Fourier trans- 
form of a Heaviside function at non-zero frequency is pro- 
portional to w"^, and the two Coulomb fields interfere 
coherently at low frequencies, explaining the frequency 
dependence of the electric field. Naturally the ZHS al- 
gorithm does not reproduce this behavior which is not 
associated to radiation. The growth of the field at low 
frequencies is an artifact of not accounting for screening 
of the field by the atoms in the medium. 

2. Behavior of the field wtth distance 

In Fig. [5] the dependence of the Fourier components 
of the field modulus with distance is shown for several 
frequencies. At sufiiciently large distances to the track 
the electric field behaves as for all frequencies. This 
is the radiation zone, the field is expected to behave as 
1/R as explained in conventional radiation theory If 
the observer is placed at very small distances compared to 
the length of the track, the situation resembles that of an 
infinite track. In this case the discussion in Section Hill 
applies. The field behaves as 1/VR at distances much 
smaller than the length of the track provided the fre- 
quency is high enough to satisfy \up\ 1 - see Eq. ((55t . 
This can be seen in the curve of the Fourier component 
ai V = \ GHz for R ^ 0.5 m. At small distances and 
sufficiently low frequencies, when \up\ <^ 1, the field be- 
comes proportional to and independent of frequency 
- see Eq. (|29p . This feature can be appreciated in Fig. [5] 
at distances below 0.1 m for the calculations at 10 and 
100 MHz. The ZHS algorithm reproduces the calculation 
provided the kR ^ 1 condition is satisfied as could be 
expected, reproducing both the 1/R and 1/VR behav- 
iors. 



B. Fourier components of the electric field in 
electromagnetic showers 

It is possible to test the ZHS algorithm in a more re- 
alistic situation. In this section we compare the Fourier 
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components of the electric field predicted by the ZHS al- 
gorithm with those obtained using the exact calculation 
in a full simulation of electromagnetic showers. 

For this purpose we have applied the exact solutions of 
the field of a track given in Eqs. (HU) and ([12]) in the ZHS 
Monte Carlo code [l^ for the simulation of electron and 
photon-induced showers in ice. The ZHS Monte Carlo 
calculates the start and end points of small sub-tracks 
of all charged particles (electrons and positrons) in an 
electromagnetic shower down to a kinetic energy thresh- 
old of ^ 100 keV. With these we can calculate the exact 
electric field produced by each single sub-track and add 
the fields up accounting for interference between differ- 
ent tracks. Since Eqs. ([TT|) and ([T^ are only valid for 
a charged particle travelling along the z axis (parallel to 
the shower axis), we perform the necessary rotations of 
Eqs. ([TT]) and (IT^ to obtain the field for a particle track 
moving along an arbitrary direction. 

Simultaneously with the exact calculation, we also ob- 
tain the field as predicted by the ZHS algorithm for ex- 
actly the same shower (i.e. the same set of tracks and 
sub-tracks). As explained above the subdivisions are such 
that the conditions in Eqs. (dH), ^ and (UHl) are fulfilled 
for all the sub-tracks in the shower. 

The result is qualitatively the same as in the case of 
single tracks. As long as the condition kR ^ 1 is fulfilled, 
the ZHS algorithm gives an accurate prediction for the 
Fourier components of the electric field with a difference 
of less than a few percent relative to those obtained with 
the exact calculation. As can be seen in Fig.[6]this occurs 
for distances to the shower axis as small as i? = 1 m 
and frequencies above ly ^ 10 MHz, well in the distance 
and frequency ranges relevant for experiments looking 
for par ticle shower induced radio pulses in dense media 

We stress here that the accuracy reported above refers 
to the approximation of using the ZHS formula applied 
to the standard subdivision of tracks in the ZHS code, 
instead of the exact expression for the radiation emitted 
by the same particle sub-tracks. By comparing the re- 
sults obtained in the Fraunhofer limit with the standard 
subdivision of tracks to those obtained with a much finer 
subdivison, it was determined that the accuracy of the 
ZHS code is ~ 10% at frequencies ~ 5 GHz, improving 
significantly at lower frequencies. We do not further ad- 
dress this uncertainty in this paper, nor the uncertainty 
due to the shower simulation itself. 

It is also worth remarking that in terms of computing 
time the exact calculation is roughly a factor ~ 5 slower 
than the calcuation performed with the ZHS algorithm. 

Since the ZHS algorithm can only be applied in a lim- 
ited range of frequencies, an accurate representation of 
the electric field in the time-domain cannot be obtained 
with an inverse Fourier transform of Eqs. and (P7|. 
The low frequency components which do not satisfy the 
condition kR ^ 1 are not accurately described by the 
ZHS algorithm as shown before. Also at very high fre- 
quencies the number of steps in which the tracks have 



to be divided in order to fulfiU Eqs. (HH), ^ and p9)) 
can become prohibitively large from the computational 
point of view what can compromise the calculation of 
the pulses below the 100 picosecond scale. In practice 
however these techniques, although not exact, provide 
fast and accurate calculations in the region of interest to 
UHE neutrino detection. For experiments with typical 
time resolutions of the order of 1 ns, and which are only 
sensitivity to frequencies from > 10 MHz to few GHz, the 
ZHS algorithm has been shown to give a very accurate 
representation of the Fourier components of the electric 
field (Fig.©. 

V. COMPARISON TO OTHER CALCULATIONS 

Several calculations of the field emitted in showers de- 
veloping in dense media can be found in the literature. 
In 47] the Finite Diference Time Domain method is used 
for calculating the field of a pancake-like shower with a 
Gaussian longitudinal development and Gaussian radial 
profile in the time-domain which is then transformed to 
the frequency-domain. In [46| using the saddle-point ap- 
proximation, an analytic equation for the calculation of 
the electric field of a charge distribution exhibiting a lon- 
gitudinal profile with a well-pronounced maximum is de- 
rived. The result is factorized into an integral accounting 
for the longitudinal variation of the charge and a form 
factor that accounts for the lateral spread of the shower, 
a procedure revisited in ^41i] for realistic showers. Assum- 
ing a Gaussian longitudinal and lateral development for 
the charge distribution both results were directly com- 
pared and turned out to be in good overall agreement as 
shown in [47|. Minor differences could be attributed to 
the form factors used. 

With the exact calculation of the electric field per- 
formed in this work, the field due to a Gaussian profile 
can also be obtained and compared to the calculations 
mentioned before. The electric current for a shower with 
a Gaussian profile is given by Eq. ([1]) with the following 
replacements: 

Z{t) = 1 (40) 

and 

P(x, t) - ^ e-'-^'W)/2.l e-'V^-'fsiz ~ vi) (41) 

The shower develops in the longitudinal direction parallel 
to the z' coordinate (shower axis), and radially along the 
x' and y' coordinates. A'^ is a normalization constant, cr; 
characterizes the width of the shower along shower axis 
and Or the corresponding lateral width. After substitut- 
ing this current in Eq. ([7]), and following the same steps 
as in Section III A[ the expression for the field is the same 
as in Eqs. (ITT|) and with the following change: 

/dt' ^ / di' dx' dy'-^ e-(-"+^")/2'^^ e-^*"/^'^? 
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FIG. 6. Fourier components of the electric field modulus 
as obtained in Monte Carlo simulations of a 10 TeV electron- 
induced shower in ice, with the exact calculation (lines) and 
with the ZHS algorithm (symbols). The field is shown for 
observers at distances (from top to bottom panels) R = 100, 
10, 1 m, placed at different observation angles with respect 
to shower maximum. In the middle panel corresponding to 
7? = 10 m, we also show the relative difference (in %) between 
the electric field modulus as obtained with the exact solution 
and with the ZHS algorithm for the various observation angles 
depicted. 



Eqs. pT|) and after the changes in Eq. (|42)) can be 
solved numerically. The ZHS algorithm can also be ap- 
plied to this situation as long as the condition fcr ^ 1 
for all distances r to the shower. The procedure consists 
on slicing the volume occupied by the bulk of the shovifcr 
in small cubes and approximating each as a track with 
constant charge given by the Gaussian distributions in 
Eq. gl]). 
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FIG. 7. Fourier components of the modulus of the electric 
field for a Gaussian charge profile given in Eq. (|4ip with ai = 
20 m and CTr = 1 m for an observer at i? = 300 m with 
respect to the peak of the Gaussian longitudinal profile. The 
observation angles are, from top to bottom, 9c, 9c + 5° , 6c + 
10°, and 6c + 20°. Fields are calculated with the saddle- 
point approach and the exact formula Eqs. ([lip and p2|) . 
The result obtained with the ZHS algorithm is on top of the 
exact calculation and it is not plotted in the Fig. for clarity. 



Comparison of the result of the exact calculation in 
this work (or the ZHS algorithm) with the saddle-point 
calculation in (46j requires knowing the form factor F for 
a Gaussian profile. F is defined in [4^ as: 



dx'dj/'ds' e-'i-^7(s',a;',?/') (43) 



with q = (w/w, kp/R) with p — {x,y) the radial position 
of the observer. Also, x' = {s' ,x',y') with s' = z' — vt' . 
R is the distance from the maximum of the shower to the 
observer. The function / represents a normalized charge 
density of the travelling pancake. Assuming a Gaussian 
for / of the form: 



1 



2'Ka': 



As') e 



(44) 



and substituting / into Eq. (|43l) . the form factor reads, 
Fiq) = e^i^'^i"-^' (45) 

Setting iV = 1, CTi = 20 m, = 1 m and i? = 300 
m with the refractive index of ice n = 1.78, the electric 
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fields for a Gaussian charge profile were calculated us- 
ing the saddle-point approach as described in [i^, and 
compared to the exact formula given in this work. The 
results are shown in Fig. [T] Since R is large the con- 
dition kR ^ 1 is satisfied for frequencies above v ^ 1 
MHz (see Eq. [36l) and the agreement between the ZHS 
and the exact calculation presented in this work is very 
good (not shown in Fig. [7] for clarity). The saddle-point 
approach is also in very good agreement with both the 
exact and ZHS calculations. The results of the field for a 
Gaussian charge profile are also in good agreement with 
those obtained in [47| using the Finite Difference Time 
Domain method. 



VI. SUMMARY AND CONCLUSIONS 

We have obtained the results of an exact calculation 
of the electric field produced by a charged particle that 
accelerates instantaneously moves at constant speed and 
intantaneously deccelerates. The results are used to ob- 
tain the approximate expression used in the ZHS algo- 
rithm to calculate the radio emission from showers in 
dense media using shower simulations. This allows the 
precise determination of the conditions necessary for this 
approximation to be valid, namely, the observer must be 
in the far field zone kR ^ 1 and the Frauhofer approxi- 
mation must apply, kL^ sit? S/i? ^ 1 for each track. The 
exact and simulated results are compared for tracks in a 
variety of circumstances to illustrate both the behavior 
of the fields and the validity of the approximation made 
in different frequency and distance ranges. 

We have shown that the range of validity of the expres- 
sions can be greatly enlarged by subdividing long charged 
particle tracks in smaller sub-tracks and adding up the 
conributions, as it is done in the ZHS algorithm. Never- 
theless the far field condition kR ^ 1 has been shown to 
be an intrinsic limit of the ZHS algorithm. The ZHS al- 
gorithm has been tested for long tracks in regions where 
the Fraunhofer approximation does not hold comparing 
it to the exact solution. The results clearly indicate that 
the ZHS algorithm reproduces the exact behavior pro- 
vided the far field condition is satisfied and the length of 
the sub-tracks is small enough for the Fraunhofer approx- 
imation to be valid. By comparing the emission with the 
exact emission from an infinite track it is shown that the 
emission calculated with the ZHS algorithm does indeed 
contain what is conventionally described as Cherenkov 
radiation. The precision of the ZHS algorithm is shown 
to be below the 2% level provided that kR > 3.7, corre- 
sponding to i? > 10 m and v > 10 MHz. By enforcing 
kR > 37 the precision improves to better than 0.01%, 
what corresponds to i? > 10 m and ly > 100 MHz, which 
are conditions met in most experimental arrangements 



trying to detect radio pulses form neutrinos interacting 
in dense media. 

The ZHS algorithm is tested for completitude when 
applied to a shower simulation. This is done compar- 
ing the result of the ZHS algorithm to that obtained 
when the ZHS expression is replaced by the exact cal- 
culation for every charged particle sub-track for exactly 
the same shower. The results indeed confirm that the 
accuracies reported for individual sub-tracks are approx- 
imately mantained in the final ZHS result. 

Finally, in order to compare to alternative calculations 
the results are compared to the saddle point approxima- 
tion. This approximation has been used to test solutions 
using the method of Finite Differences to solve Maxwell's 
Equations directly in the Time Domain (FDTD) using a 
simplified shower front based on gaussian distributions 
in the shower plane and in time. The comparison of the 
saddle point approximation to both the exact solution 
and the ZHS algorithm give compatible results confirm- 
ing that both approaches reproduce the radiation emitted 
from showers. 

The results presented here in summary confirm that 
the ZHS algorithm can be used to describe most practi- 
cal applications to detect pulses emitted from high en- 
ergy showers produced in dense media by neutrinos. The 
approach only begins to show significant discrepancies 
when the observer is at distances comparable to the lat- 
eral dimensions of the shower (< 1 m in ice). Since the 
typical distance between antennas in experiments such 
as the Askaryan Radio Array (ARA) [22] is 10 - 100 
m, we expect the results to be accurate enough in most 
practical situations. 

The results obtained are also of interest for the ap- 
plication of the ZHS algorithm to calculate pulses from 
EAS. Indeed the application of the method of track sub- 
division has been applied in ZHAireS [i^, a recent code 
developed for this purpose, by making track subdivisions 
that are forced to satisfy the Fraunhofer condition. 
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